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Abstract. Based on the Chetaev theorem on stable dynamical trajectories in the presence 
of dissipative forces, we obtain the generalized condition for stability of relativistic classical 
Hamiltonian systems (with an invariant evolution parameter) in the form of the Stueckelberg 
equation. As is known, this equation is the basis of a competing paradigm known as 
parametrized relativistic quantum mechanics (pRQM). 

It is shown that the energy of dissipative forces, which generate the Chetaev generalized 
condition of stability, coincides exactly with Bohmian relativistic quantum potential. Within 
the framework of Bohmian RQM supplemented by the generalized Chetaev theorem and on the 
basis of the principle of least action for dissipative forces, we show that the squared amplitude 
of a wave function in the Stueckelberg equation is equivalent to the probability density function 
for the number of particle trajectories, relative to which the velocity and the position of the 
particle are not hidden parameters. 

The conditions for reasonableness of trajectory interpretation of pRQM are discussed. Based 
on analysis of a general formalism for vacuum-flavor mixing of neutrino within the context of 
the standard and pRQM models we show that the corresponding expressions for the probability 
of transition from one neutrino flavor to another differ appreciably, but they are experimentally 
testable: the estimations of absolute value for neutrino mass based on modern experimental 
data for solar and atmospheric neutrinos show that the pRQM results have a preference. It 
is noted that the selection criterion of mass solutions relies on proximity between the average 
size of condensed neutrino clouds, which is described by the Muraki formula (29th ICRC, 2005) 
and depends on the neutrino mass, and the average size of typical observed void structure (dark 
matter -I- hydrogen gas), which plays the role of characteristic dimension of large-scale structure 
of the Universe. 



1. Introduction 

In this paper we continue to consider the question, which is formulated in [1, 2] in the following 
rather strict and paradoxical form: "Are the so-called quantization conditions that are imposed 
on the corresponding spectrum of a dynamical system possible in principle in classical mechanics, 
analogously to what is taking place in quantum mechanics?" As is shown in [1, 2], the answer to 
this question is positive and has been given more than 70 years ago by the Russian mathematician 
and mechanician N.G. Chetaev in his article "On stable trajectories in dynamics" [3, 4, 5, 6]. 
What is more, based on the Chetaev theorem on stable dynamical trajectories in the presence 



of dissipative forces [3, 4, 5, 6], we obtained the generalized condition for stability of Hamilton 
systems in the form of the Schrodinger equation [1, 2]. 

At the same time, it is necessary to remind that the essential difference of the Chetaev 
quantum mechanics with trajectories from the known alternative quantum theories [7-16] is the 
unconditional identity of the probability density function of the number of particle trajectories 
(obtained by the principle of least action of perturbations (see (28) in [1, 2]) and the probability 
density function for a particle to be in a certain place of the configuration space (obtained by 
the Bohm-Madelung equation of continuity (see (25) in [1, 2])). This exclusively important fact 
emphasizes naturally the physical identity of the probabilistic and trajectory interpretations of 
nonrelativistic quantum mechanics. 

In this paper we show that similar results can be obtained in the framework of relativistic 
classical mechanics. In other words, based on generalized stability conditions for stability of 
Hamilton relativistic systems (with invariant evolution parameter) we obtained the relativistic 
Schrodinger equation or, more exactly, the so-called Stueckelberg equation, which as is known 
is the basis of a competing paradigm known as parametrized relativistic quantum mechanics 



On the other hand, it is obvious that the fact of appearance of an alternative quantum 
theory causes immediately a number of questions connected, in the first place, with comparison 
between results of alternative and standard quantum theories obtained in the same experiments. 
As is exactly mentioned by Adler and Bassi [18], "...but apart from this history, there is 
another important motivation for considering modifications of quantum theory. This is to 
give a quantitative meaning to experiments testing quantum theory, by having an alternative 
theory, making predictions that differ from those of standard quantum theory, to which these 
experiments can be compared" . 

Thus, the main purpose of the present work is, on the one hand, derivation of the 
Schrodinger relativistic equation or, more exactly, so-called the Stueckelberg equation based on 
the generalized condition for stability of Hamilton relativistic systems (with invariant evolution 
parameter) and, on the other hand, quantitative comparison between results of alternative 
(pRQM) and standard quantum theories by analysis of general formalism for vacuum-flavor 
mixing of neutrino. 

2. The emergent Stueckelberg equation in relativistic classical mechanics 

In what follows, we present a generalization of the Chetaev theorem on stable trajectories in 
dynamics to the case of relativistic Hamiltonian K{s, q^,Pi), where s is the evolution parameter 
of a system i = 0, 1,2,3. With this purpose, we consider a material system, in which and 
Pi are the generalized 4-space coordinates and momenta of the holonomic system in the field of 
potential forces which admit a force function of the form U{q^). 

In the general case where the action S is the explicit function of the evolution parameter s, 
the total integral of the Hamilton-Jacobi differential equation corresponding to a given system 
takes the form: 



where /3* are the new constants of integration. The possible motions of the mechanical system 
are defined by various values of the constants a* and 



(pRQM) [17]. 




= dS/dui, Pi = dS/dq\ i = 0, 3 




We call a motion of the material system, whose stability will be studied, the unperturbed 
motion. First, we consider the stability of such motion with respect to the variables under a 
perturbation of only the initial values of these variables (i.e. of values of the constants a* and /3*) 
without disturbing forces. Omitting the details of the derivation which are given in [3, 4, 6, 19], 
we present the necessary Chetaev stability condition in this case as follows: 

A = —dicfS = 0, (3) 
m 

We now complicate the problem. Let a really moving material system undergo the action 
of both the forces with force function U, theoretically considered above, and the unknown 
perturbation (dissipative) forces which are assumed to be potential and admit the dissipative 
force function Q. Then the real motion of the material system occurs in the field of forces with 
general force function U* = U -\-Q. In this case, the real motion of the system does not coincide 
obviously at all with the theoretical one (without perturbation). 

If we conserve the statement of the problem on the stability of real unperturbed motions in 
the theoretical field of forces with the function [/ at a perturbation of only initial data as above, 
then the necessary requirement of stability in the first approximation, e.g., in form (3), will not 
be efficient in the general case, because the new function S is unknown (as well as Q). However, 
it turns out that we can determine such conditions of stability, which are implicitly independent 
of the form of the unknown action functions S and the potential Q. Thus, we are based on the 
requirement of stability of form (3), by assuming the conditions of its existence (correctness, 
etc.) for real motions to be satisfied. In relation (3), we now replace the function by a new 
function A defined by the equality: 

V' = ^exp(iA:S), (4) 

where A; is a constant; ^ is a function of the generalized coordinates q and the evolution parameter 
s. 

The introduction of a real wave function, like the de Broglie "pilot-wave" [20], is extremely 
necessary from the physical viewpoint because of the following non-trivial reason. Since the 
dynamics of a physical system must undoubtedly conserve the Hamilton form of the equations 
of motion, the main "task" of such a real wave consists in the exact compensation of the action 
of dissipative forces, which are generated by the perturbation energy Q. We will show below 
that, in this case, such a procedure makes it possible not only to conserve the Hamilton form 
of the dynamics of a physical system, but allows one to determine the character of an analytic 
dependence of the energy Q of disturbing forces on the wave function amplitude A. Then relation 
(4) yields: 

and, hence, relation (3) looks like 

= 0. (6) 

On the other hand, we can write the Hamilton-Jacobi equations for a perturbed motion in 
the general case where the Hamiltonian H depends explicitly on the evolution parameter s: 

^ + —diSd's + ^7 + g = 0, 

OS Zm 

or with the help of (4-6) 
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where dS/ds can be determined with the help of (4). Adding relations (6) and (7), we obtain 
the necessary condition of stability (in the first approximation) in the form: 

^ did"^ - i^did'A l^d^A ( ^dii^ - \diA 
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ikAi/j 



(8) 

[Ads^p - i^dsA] -U-Q = 0. 



In this place, we need to use a procedure for the compensation of the action of dissipative 
forces, which are generated by the perturbation energy Q in order to conserve the Hamilton 
form of the dynamics of a physical system (8). It is obvious that equality (8) will not contain 
Q, if the amplitude A is determined from the equation: 

^ did'A + -^d,Ad's-^dsA + Q = o, (9) 



2k^A ' kA ' ikA 

which splits into two equations 

dsA = --diAd'S, (11) 
m 

after the separation of the real and imaginary parts. Here Q is the dissipation energy. Thus, if 
the properties of disturbing forces satisfy conditions (10) and (11), then the necessary condition 
of stability (8) has the form of a differential equation of the " Stueckelberg" type: 

^a.v = -^W + c/v, (12) 

where g'^(s) gives the position of the physical system, whose possible trajectories in the 4- 
dimensional configuration space q = {qo,,Q3) are a solution of the system of the so-called 
governing equations [21]: 

Vi = -djS, (13) 
m 

where S is the phase of the wave function. 

In other words, we obtained the following result: Eq. (3), which corresponds to the Chetaev 
stability condition, is transformed into an equation of the "Stueckelberg" type (12) with the 
use of transformation (4). It should be noted that in spite of similarity of the Stueckelberg and 
Schrodinger equations they are different. Evolution parameter differs from the time [17, 22]. It 
is obvious that, in the class of equations of type (12), the single-valued, finite, and continuous 
solutions for the function tp in the stationary case are admissible only for the eigenvalues of the 
total energy E. Hence, the given stability of real motions takes place only for these values of 
the total energy E. 

We now present a short analysis of the obtained results. It is known that, according to one of 
the theorems of stability theory [6] , only two types of forces - dissipative and gyroscopic do not 
break the stability (if it is present) of a non-disturbed motion of holonomic mechanical systems. 
Therefore, by introducing a dissipative perturbation Q into the Hamilton-Jacobi Eq. (7) and by 



taking simultaneously the stability condition (3) into account, we reasonably expect that, under 
condition of the conservation of the stability of a mechanical system, it is possible to get a real 
functional dependence of the dissipation energy Q on characteristics of the wave function tJj (4). 
Indeed, having obtained the condition for the stability of trajectories of a dynamical system in 
the form of an equation of the " Stueckelberg" type (12), we established not only the physical 
sense of the perturbation energy Q, but showed also that it is a function of the amplitude of the 
wave function ip (4) and takes form (10). It is a very important result. The subsequent content 
of the article is, as will be clear in what follows, a direct consequence of this result. 

It is easy to show that namely the conclusion about the dissipative nature together with the 
simultaneous determination of a functional dependence of the energy of disturbing forces on the 
form (but not on the magnitude ) of the amplitude of the wave function -0 allow one to generalize 
an equation of the "Stueckelberg" type (12) to the case where the condition of stability (3) is 
not fulfilled, i.e., A 7^ 0. Let us show this. 

Obvious analysis of Eqs. (3) and (10)-(13) shows that with allowance for dimensions the 
expression: 

is, in general case, the variations of particle kinetic energy predetermined accordingly by 
variations of its momentum. 

Now, if to add and simultaneously to subtract the complex expression (ie) and to substitute 
also Eq. (4) in the left-hand side of the Hamilton- Jacobi equation (7), we get the generalized 
equation corresponding to the extended Eq. (8): 



2mk'^ijj 2mk'^A 2mk mk'^A A y ^^^^ 

Repeating the ideology of derivation of Eq. (12) it is easy from Eq .(15) to obtain an equation 
of the "Stueckelberg" type, which is formally identical to Eq. (12), but already under a more 
general condition imposed on the disturbing energy and the wave function: 

8,A = -Raid's - diAB'S. (17) 

Physical sense of (17) for the wave function amplitude consists in the fact that with regard for 
the formulas for a classical velocity Vi = diS/m and the probability density P{q,t) = [A{q,t)]'^ 
(the substantiation of this formula will be given below) in the configuration space it can be 
easily transformed into the equation of continuity, which represents the invariability of the total 
number of "particles" (phase points) or, in other words, the probability conservation law. 

We now return to our problem of quantization on the basis of the simplest example. Let 
us consider a material point with mass m in the field of conservative forces with the force 
function U, which depends, in the general case, on the evolution parameter. The problem on 
the stability of motions of such a point will be posed in the Cartesian coordinate system qq, qi, 
q2, qs- Denoting the momenta along the axes by po, pi, P2, Ps, respectively, we obtain known 
expression for the Hamiltonian: : 



K 



2m 



PiP ■ 



(18) 



In this case, conditions (16) and (17) for the structure of disturbing and compensating forces 
admit the relations: 



Q 



did' A 
2m A 



k = l/h , 



A 
2m 



V — 



E 



dA p' 



dqf ^ dqi m 

and the differential equation (12) defining stable motions takes the form: 
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(19) 
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ds 2m 

which coincides with the Stueckelberg equation. In our case, this equation restricts the choice 
of the integration constants in the full Hamilton- Jacobi integral. In what follows, we call Eq. (21) 
the Stueckelberg-Chetaev equation, by emphasizing the specific feature of its origin. 

It is of interest to consider the case related to the inverse substitution of the wave function 
(4) in the Schrodinger equation (21), that generates an equivalent system of equations known 
as the Bohm-Madelung system of equations [7, 8, 9]: 
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It is very important that the last term in Eq. (23), which is the "quantum" potential of the 
so-called Bohm -(/'-field [7, 8, 9, 23] in the interpretation of Ref. [7], coincides exactly with the 
dissipation energy Q in (19). At the same time, Eq. (22) is identical to the condition for dA/ds 
in (17) and (20). 

If we make substitution of the type: 



P{q,s)=ipr = [A{q,s)f, 
then Eqs. (22) and (23) can be rewritten in the form: 
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(24) 



(25) 



(26) 



Here, Eq. (25) has a clear physical sense: P(q, s) is the probability density to find a particle 
in a certain place of the 4-space, and diS/m is, according to (13), the 4- velocity of this particle. 
In other words, Eq. (25) is nothing but the equation of continuity which indicates how the 
probability density P{q, s) moves according to the laws of classical mechanics with a classical 
velocity at every point. 

On the other hand, we can show that P{q, s) is also the probability density function for the 
number of particle trajectories, that is substantiated in the following way. We assume that the 
influence of disturbing forces generated by the potential Q on a wave packet at an arbitrary point 



of the configuration space is proportional to the density of particle trajectories (iptp* = A^) at 
this point. This implies that the disturbing forces do not practically perturb the packet, if the 
relation 



QiPi;*dV min, / V?/'*^^ = 1, (27) 



is satisfied, where dV stands for an element of the configuration space volume. This means, 
in its turn, that the disturbing forces admit the absolute stability on the whole set of motions 
in the configuration space only if condition (27) is satisfied or, in other words, if the following 
obvious condition of the equivalent variational problem (for Q) fulfills: 

5 j Qi)^*dV = = 0. (28) 

The variational principle (28) is, in essence, the principle of least action of a 
perturbation. Below, we call it the principle of least action of perturbations by Chetaev [1]. 
Using the previous notations, we can write the following equality for Q: 

BS 1 

Q = -dsS-U-T = -—-U--diSd'S. (29) 
On the other hand, if (4) holds true, it is easy to show that: 

Then it is necessary to perform the following subsequent substitutions. First, we substitute 
relation (30) in (29) and then introduce the result in the equation corresponding to the variational 
principle (28). 

As a result of the indicated procedure of substitutions, we obtain the relation which is exactly 
equal to Eq. (15), i.e., to the extended equation (with regard for (3), (14), (15)] of type (8). 
Hence, the structural expressions and the necessary condition of stability which follow from it 
coincide with (16), (17), and (12), respectively. This means that, on the basis of the Chetaev 
variational principle (28), we get the independent confirmation of the fact that the physical 
nature of P(q, t) reflects really not only the traditional notion of the probability density for 
a particle to be at a certain place of the space according to the Bohm-Madelung equation of 
continuity (25) but plays also the role adequate to that of the probability density of the number 
of particle trajectories. 



3. Neutrino oscillation measurements as verifiability criterion of different concepts 
of temporal evolution 

Here the natural question raises: "What of the two equations realizing ideology of standard or 
parametrized RQM is substantiated and valid within the framework of modern experimental 
relativistic dynamics or, in other words, what of the two concepts of temporal evolution reflects 
nature of physical reality of our world more sufficiently?" It is interesting, that one of variants 
of the theoretical formulation of such a problem is stated for the first time by Fanchi [17, 22]. 
Because we will use the results of this work below, let us briefiy consider the formulation of this 
problem. 

Following [22], we consider the temporal model of vacuum-flavor mixing neutrino within the 
framework of the standard and pRQM models. At first, this is caused by the fact that, according 
to [22], differences between the standard and pRQM models are highlighted by developing a 
general formalism for vacuum-flavor mixing and then applying the formalism within the context 
of each theory. Secondly, the resulting two expressions for the probability of transition from 



one neutrino flavor to another are experimentally testable that is the one of purposes of our 
research, which is given in next sections. 

Considering the temporal model of vacuum-flavor mixing, we imply the case of two-state 
mixing, i.e., we confine ourself to consideration of the two mass states {| Uj)} and two neutrino 
flavor states {| fa)}, which may be written as the two-component column vectors: 



{K)} 



(31) 



The temporal evolution equation for state may be written in terms of the temporal evolution 
operator T as: 



with the formal solution 
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It is obvious [22], that the eigenvalue Tj of the temporal evolution operator T depends on 
the context of theory. In the standard theory, T-,- is energy Ej of state j, while in pRQM it is 
the eigenvalue Kj of the mass operator for mass state j. Explicit free-particle expressions for 
the eigenvalues are: 



(T 



j ) standard 



,1/2 



'^'Standard 



(34) 



and 



(Tj)pPQM = Kj - 



(35) 



2m i 



(ujj/c^) - kj ■ kj , {t)prqm 



Here 



Wi, k 



I denotes the energy-momentum four-vector in laboratory frame and 



mi- 



ls the mass of state for a metric with signature {+, ■ 



-}. According [22], in pRQM, 



all four component of the energy-momentum four-vector are observable and nij is function of 
statistical values of the energy-momentum four-vector. In this case the energy u)j of state j is 
treated in pRQM as an independent variable and does not satisfy an equation like (34), except 
in a statistical sense. By contrast rrij is an input parameter in the standard model. In this 
case the standard picture of neutrino oscillation by flavor mixing in vacuum assumes that the 
three- momentum is constant for all states and has magnitude = h?k ■ k. 

It is known that a neutrino oscillation process is descried by hypothesizing that the mass 
basis I Uj) is related to the flavor basis | z^o) by a unitary transformation U such that 



{Wa)} = U{W,)} 



cos t> 
— sin( 



smf 
cos ( 



(36) 



and 6 is the mixing angle of mass states in vacuum. For two flavors of neutrinos and the unitary 
transformation U we have the r-dependent state: 



(37) 



The probability of transition from I'eiT) to is the square of the magnitude of the transition 
matrix element 



Thus, the survival probability P(fe — ?• i^e) has the form 

P {ue -^Ue) = l- P{ue ^v^,) = l- sin^ 26 sin^ (39) 

It is obvious that, that standard model and pRQM model according to (34-35) differ in the 
explicit forms of Tj and r. Therefore, the expression for the survival probability P{ye ~^ ^e)i 
obtained in the framework of standard model and pRQM model, will have the next form: 
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(40) 



where the coordinate time t is replaced in standard model by the distance L from the neutrino 
source to the detector, so that t = L/c. Note that for ultra-relativistic neutrino it is assumed 
that that pi ~ P2 = P-, and that p/c » 1712, mi, and p ^ Ey/c, where Ey = mic^/(l — /3^)^/^ is 
the energy of the incident neutrino (/3 = v/c): 
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(41) 



Based on well-known experiments on solar and reactor neutrino detection [24] one can to 
conclude that existence of neutrino oscillations is indisputably proved and to assume without 
loss of generality that the mass difference between neutrino flavor states is very small, i.e., 
mi ~ m2. Thus, it is easy to see that in the case of ultra-relativistic neutrino (i.e., /3 — )• 1) the 
dynamical factors o-standard and apRQM are connected by the following simple expression: 



OlpRQM = 2a Standard- 
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where Am^i = m^ — mf. 

In other words, the ratio of dynamical factors does not equal to 1 because of the difference 
in operators Kj between standard paradigm and pRQM. It was pointed out in Fanchi [22] 
that the pRQM mass operator is the direct result of the probabilistic basis of the theory. 
Consequently according to [22] , the search for massive neutrinos is a test not only of the flavor- 
mixing hypothesis, but also the validity of standard model and pRQM. 



4. Neutrino mass matrix via see-saw mechanism 

It is obvious, that to analyse a general formalism for vacuum- flavor mixing of neutrino within 
the context of the standard and pRQM models we need the estimations of absolute value for 
neutrino mass based on modern experimental data for solar and atmospheric neutrinos. With 
that end in view we consider below the problem of determining the neutrino masses by use of 
special method for construction of neutrino mass matrix, which was described for the first time 
in [25], and advantages of the mass-squared differences and mixing parameters obtained from 
neutrino oscillation experiments. 

It is known that one of the interesting mechanisms which can generate a small neutrino mass 
is the see-saw mechanism, in which the right-handed neutrino vpt has a large Majorana mass 
Mat and the left-handed neutrino vl obtain a mass through leakage of the order of (m/Mjv) 
with ni is the Dirac mass. According to see-saw mechanism [26], the neutrino mass matrix M^j 
is given by: 

« -MdM^^mI, 

where Mr, and are the Dirac and Majorana mass matrix respectively. 

In order to obtain neutrino mass matrix Mi, that can give correct predictions on mass-squared 
differences and mixing parameters, Damanik [25] modified the the neutrino mass matrix Mi, by 
introducing one parameter 6 to perturb the diagonal elements of My. In the way that the 
perturbed mass matrix satisfies the requirement Tr (M^) = 3P. At the same time, the neutrino 
mass matrix M^ can be written in the following form [25]: 

/ P + 25 Q Q \ 
My=\ Q P-6 Q . (44) 
\ Q Q P-6 J 

The eigenvalues of the neutrino mass matrix M^ in (44) are: 



,,,^P^<l,l^VMEM±W, ,45, 

P3 = P-Q-S. (46) 
If the neutrino mass matrices M^ in (44) is diagonalized by mixing matrix V in the equation 
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then we obtain 



and neutrino mass as follow: 
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(51) 



m3 = P - Q - 6. (52) 

One can see that the obtained neutrino masses in Damanik's scenario [25] is an inverted 
hierarchy with masses: | ms |<| mi \<\ m2 |. 

To solve the system of the equations (49)- (52) with respect to the values of neutrino masses 
it is necessary to supplement this system with the known values of mass-squared differences 
and mixing parameters, which obtained from neutrino oscillation experiments. From the recent 
global analysis of the neutrino oscillation data we refer to the following best-fit values of solar 
neutrino mass-squared differences [27]: 

Amli = (r.btoil) X 10"^ eV"^, (53) 

with 

tan^^ia = 0.452l°:0|, (54) 
and for the atmospheric neutrino mass-squared differences [28]: 

Am|2 = (2.32l°;J|) x 10^^ eV"^. (55) 

If 9 is the 012 in (54), we have from (49) we have 6 = —0.050909 Q. If we insert this S value 
into (50-52) then we obtain the neutrino mass as follow: 

mi = P- 1.052610 Q, (56) 
m2 = P 2.001700 Q, (57) 

= P - 0.949910 Q. (58) 

Using the advantages of experimental data of neutrino oscillations in Eqs. (53) and (55) from 
Eqs.(56,57), we have the neutrino masses Table 1) obtained within the framework of standard 
model and pRQM model. 

Table 1. Neutrino mass |mj| as function of dynamical factors a standard and otpRQM- 
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0.172499 



Thus, we estimated the neutrino masses within the framework of both models. So long as 
these estimations satisfy the known restriction on the total neutrino mass: 

^ ma < 0.6eV, Va = Ve, Vr, (59) 

it is impossible to answer question, what of the two representations realizing the ideology 
of standard or parametrized RQM reflects the nature of physical reality of our world more 
sufficiently. A valid answer to this question can be given only when the experimental values of 
the neutrino masses will be obtained. 

Below we will attempt to solve this problem even if at qualitative level using modern 
cosmological observations. 



5. The dimension of the neutrino cloud and large-scale structure of Universe 

In this section we compare the theoretical value of void structure of Universe obtained by use 
of foregoing alternative values of neutrino masses with the experimental value of average size 
of large-scale structure of Universe obtained on the basis of analysis of so-called the matter 
power spectrum of density fluctuation in space-time. The main purpose of such a comparison 
is the reveal of dominant role one of the two physical concepts of temporal evolution, which are 
considered above. 

So, according to Muraki [30, 31], the possible role of massive neutrinos in our Universe is 
as follows. Massive neutrinos are attracted by the gravitational force. However, because of the 
degeneracy force, they will not collapse to form a supper black hall. The Jeans mass of neutrinos 
grows and fluctuation of the baryonic matter will be erased. The fluctuation of the baryonic 
matter can grow only at the border between neutrino bubbles. Consequently, it is possible to 
assume that the observed void structure of the Universe was formed by the neutrino degenerate 
force and grows of the density fluctuation at the border between voids. 

Using the analogy with neutron star, whose radius ^ can be obtained by the balance between 
the degenerate force of neutrons and the gravitational force, Muraki [30, 31] in bounds of his 
script (scenario) got the next expression for the radius of a neutrino cloud: 



where the neutrino cloud mass is 



= (4/3) Txrl X 110 X 10^ [m?] x [eV] x 1.7 x IQ-^^ [kg] . (61) 

However, in this scenario it is necessary to take into account also the effect of baryonic 
matter and dark matter [30, 31]. According to the last analysis of the WMAP7 data, total 
matter fraction 17^/ is = 26.6 it 2.9% [32], while the energy density of all types of neutrinos 
is rjjy < 0.014 [33] ^. That means that the mass of (dark matter + hydrogen gas) is 19 times 
heavier than the mass of neutrinos. Therefore if we replace llOniu by 2090mi/ and at the same 
time take into account the contribution to the gravitational mass from other neutrinos, we will 
obtain modern expression of Muraki relation for the neutrino cloud diameter: 

d^ = 7.1 [Mpc] X [gy] _ (g2) 

^ The canonical radius of a neutron star [29] can be obtained by the balance between the degenerate force 
of neutrons and the gravitational force [30, 31]: r„ = 1.2(ft^/G)(l/M„)l/3(l/m„)**''^ ~ 10 km, where M„ 
corresponds to a 1.4 solar mass. 

^ In terms of neutrino mass this restriction < 0.014) means that "Y^mv < 0.6 eV [24, 33], which is satisfied 
by the values of triplet neutrino masses listed in Table 1. 
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Figure 1. The matter power spectrum P{k) predicted by the best-fit ACDM cosmological 
model (continuous curve) and how neutrino masses affect it (dashed curves). Measurements 
at different scales have been performed with different techniques, that slightly overlap. The 
data points do not show the overall uncertainty that plagues galaxy surveys (SDSS, 2dF) 
at intermediate scales and especially Lyman-a data at smaller scales i.e. at larger k. Here 
fu = Pu/pdm, where pi, - neutrino density, pdm ~ density of dark matter. The wide vertical 
stripe identifies an interval maximum best-fit ACDM cosmological model (continuous curve), 
corresponding to the range of wavelength values A = 20 500 [/i/Mpc]. Adopted from [24]. 



Using this expression it is easy to estimate the average diameter of neutrino cloud for different 
neutrino masses predetermined by the corresponding dynamical factors astandard and apRQu 
(see Table 1): 



m 



pRQM 



m 



V 



Standard 



= 0.185461ey ^ d^J^'^A/ ^ goMpc, (63) 
0.131141 eV d^t^'^'^^rd ^ i50Mpc. (64) 



On the other hand, it is known, that the average size of typical observed void structure 
(dark matter -|- hydrogen gas), which plays the role of characteristic dimension of large-scale 
structure of the Universe, are predetermined by so-called of the matter power spectrum of density 
fluctuation in space-time (Fig. 1). 

It is easy to show (see Figure 1 [24]), that the range of large-scale structure (LSS) set by the 
range of the wavelength A = 20 500 [h~^Mpc\ {h = 0.704 [32]) is described by the average 
characteristic dimension of LSS of the Universe: 

(Al55)~90 Mpc. (65) 

This is confirmed also by the known results of Tegmark et al. [34] . 

Comparative analysis of expressions (65) and (63,64) shows that the characteristic dimensions 
of LSS (65) and neutrino cloud (63) practically coincide: 



{Xlss) ~ 



(66) 



This result shows that the physical concept of temporal evolution described by Stueckelberg 
equation is more adequate. 

6. Conclusion 

Certainly, it is impossible to consider obtained result as a proof, but it is a strong argument in 
favour of the Stueckelberg concept of temporal evolution, which reflects the nature of physical 
reality of our world more sufficiently than the standard (Dirac) concept. It is clear that direct 
proof in favour of one or the other physical concept of temporal evolution will be possible only 
when the experimental value of neutrino mass will be determined. 

In this sense, it is appropriate to note here that "... time and again a breakthrough in the 
field of new physics was related to neutrino. The first theory of /3-decay proposed by Fermi was 
based on the Pauli hypothesis of existence of neutrino. Phenomenological (V-A) theory of weak 
interactions has arisen from the two component neutrino theory considered by Landau, Lee, 
Yang and Salam. Evidence in favour of the Standard Model of weak interaction of Glashow, 
Weinberg and Salam had been received on a neutrino beam (discovery of neutral currents)" [35]. 
At the present time, we observe a new phenomenon - neutrino oscillations, which result from 
temporal evolution of massive neutrinos and indicate new physics beyond the Standard Model. 
Undoubtedly, that discovery of neutrino mass in the near future (for example, by research of 
high-energy edge of the tritium /3-spectrum) can become the basis for reconsideration not only of 
the key issues of quantum theory, but a holistic understanding of the nature of physical reality. 
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